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I  INTRODUCTION 

Development  of  advanced,  high-speed  aircraft,  missiles  and 
re-entry  vehicles  requires  specialized  application  of  scientific 
knowledge  from  many  fields.  One  key  technical  area  is  aero- 
ballistics  research.  The  present  paper  deals  with  a  specific 
topic  in  this  field,  namely,  effects  of  surface  roughness,  and 
roughness  plus  waviness,  on  the  surface  shear  stress  and  velocity 
profile  of  a  compressible,  turbulent  boundary  layer  in  adiabatic 
flow. 

Turbulent  boundary- layer  flows  over  smooth  surfaces  have 
been  studied  in  great  detail  and  numerous  texts  and  papers 
documenting  advances  in  the  state  of  the  art  have  been  written 
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(e.g.,  refs.  (1)  through  (7)).  Research  on  effects  of  roughness 
on  turbulent  boundary-layer  characteristics  has  not  received 
the  same  overall  level  of  attention.  However,  a  core  of 
notable  experimental  results  has  been  generated  in  this  area 

and  a  brief  review  is  outlined  below. 

(8) 

The  works  of  Nikuradse v  ',  for  incompressible  flow,  and 
(9) 

Goddard  ,  for  compressible  flow,  provided  an  early  basis 
concerning  effects  of  roughness  on  the  skin  friction  and  velocity 
profile  of  a  fully  developed  turbulent  boundary  layer.  Results 
substantiating  and  extending  the  findings  of  these  classic 
studies  have  been  reported  by  other  investigators  (e.g.,  refs. 

(10)  to  (15)).  Combined  interactions  of  surface  roughness, 
skin  friction,  and  heat  ransfer  have  also  been  studied  (e.g., 
refs.  (16)  to  (22)).  Based  on  such  experimental  evidence,  semi- 
empirical  prediction  techniques  for  turbulent  heat-transfer 
rates  in  the  presence  of  surface  roughness  have  been  formulated 
(e.g.,  refs.  (23)  to  (29)). 

Along  similar  lines,  effects  of  surface  waviness  on  turbulent 
boundary-layer  characteristics  have  been  studied  both  theoretically 
and  experimentally  (e.g.,  refs.  (30)  to  (35));  in  all  cases, 
however,  the  wavy  surfaces  under  consideration  were  aerodynamically 

smooth,  i.o„,  with  respect  to  superimposed  roughness. 

Turbulent  boundary-layer  flows  over  surfaces  which  possess 
a  roughness  pattern,  simultaneously  superimposed  on  one  or  more 
periodic  waveforms  (i.e.,  rough/wavy  walls),  have  yet  to  be 
investigated.  The  present  paper  addresses  this  problem  for 
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the  case  of  adiabatic  flow.  This  research  was  motivated  by 
an  ongoing  missile  development  program,  from  which  both  the 
need  for  such  information,  and  its  unavailability,  were 
identified. 

Current  submarine- launched  ballistic  missiles  possess 
motor  cases  made  of  a  fiber  glass-wound  material.  Surface  ridges, 
or  grooves,  plus  any  surface  undulations  created  during  the 
fabrication  process  are  everywhere  transverse  to  the  local  flow 
direction.  Such  multistage-boost  vehicles  have  a  large  per¬ 
centage  of  their  surface  area  covered  with  such  waves  and 
grooves . 

During  recent  flight  tests,  range  reductions  were  experi¬ 
enced  which  could  not  be  accounted  for  by  original  drag 
calculations.  It  was  felt  that  the  random  rough/wavy  surface 
pattern  of  these  motor  cases  was,  in  some  manner,  amplifying 
skin-friction  drag  levels  beyond  estimates  based  strictly  on 
the  smallest  roughness  dimension,  the  filament  radius. 

An  experimental  program  was  thereby  formulated  to  investi¬ 
gate  turbulent  flows  over  such  complex  surface  patterns. 

II  ROUGHNESS  MODELS 

Surface  contour  traces  taken  from  an  actual  motor  case 
section  showed  the  existence  of  three  dominant  features,  i.e., 
physical  scales,  (1)  a  roughness  scale,  (2)  a  short  wavelength 
scale,  and  (3)  a  long  wavelength  scale.  In  an  attempt  to  simulate 
these  features,  several  rough/wavy  patterns  were  created. 

Figure  1  schematically  shows  the  most  complex  simulation 
fabricated  for  testing,  wherein  roughness  was  superimposed  on 
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a  short  wavelength  waveform,  both  of  which  were  superimposed 
on  a  long  wavelength  waveform.  Table  1  summarizes  the  roughress 
models  tested. 

As  can  be  seen,  three  general  categories  of  models  were 
tested,  sand-grain  and  machined  roughnesses,  plus  an  actual 
mold  of  the  roughness  pattern  taken  directly  from  a  full-scale 
motor  casing. 

Sand-grain  models  were  fabricated  by  bonding  standard, 
uniform-grit,  sandpaper  sheets  on  to  aluminum  plates.  Machined 
models  were  actually  fabricated  by  impressing  the  desired  surface 
patterns  on  aluminum  plates;  this  technique  required  the 
machining  of  steel  rollers  with  the  desired  patterns  which,  in 
turn,  were  repeatedly  traversed  across  the  softer  aluminum 
surfaces.  Roller  widths  were  many  times  the  roughness  dimensions, 
and,  where  appropriate,  were  equal  to  several  short  wavelengths, 
or  one  long  wavelength.  Amplitude-to-wave length  ratios  generated 
by  this  technique  were  the  order  of  0,010  for  the  short  waves 
and  0.005  for  the  long  waves.  All  dimensions  shown  in  Table  1 
were  verified  through  optical/photographic  techniques  and 
surface  contour  traces  (Pipe  Machinery  Co.  contour  reader,  with 
wedge-tipped  stylus  0.001  inch  in  thickness;  system  accuracy 
quoted  at  ±0.2  mil). 

Three  categories  of  surface  roughness  were  included  in 
order  to  meet  the  following  requirements:  (1)  the  smooth  wall 
model  would  serve  as  a  reference  case,  (2)  both  the  smooth  wall 
and  sand-grain  models  would  serve  to  provide  a  check  on  present 
experimental  apparatus/instrumentation  and  testing  techniques 
tnrough  comparisons  of  shear  and  profile  data  with  previously 
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published  results,  (3)  sand-grain  and  molded  model  shear  data, 
in  combination,  would  serve  to  define  equivalent  sand-grain 
roughnesses  for  the  motor  case  material,  and  (4)  the  machined 
roughness  samples  woul  serve  to’ show  effects  of  roughness, 
with  and  without  a  single  periodic  waveform,  as  well  as  effects 
of  simultaneously  combining  roughness  with  two  periodic  wave¬ 
forms  of  differing  scales. 

Ill  EXPERIMENTAL  APPARATUS  AND  TECHNIQUES 

All  tests  were  run  in  the  Naval  Ordnance  Laboratory's  (NOL) 

Boundary  Layer  Channel^),  a  vertical,  asymmetric,  variable 

Mach  number  facility,  comprised  of  a  flat  nozzle  wall  (test 

surface)  and  a  contourable  nozzle  wall  (for  variations  in 

free-stream  Mach  number ,  pressure  gradient,  etc.).  This 

facility  has  been  used  for  many  detailed  studies  of  smooth 

wall  turbulent  boundary  layers,  e.g.,  references  (37)  to  (39), 

and  its  capabilities  are  well  documented.  For  this  program, 

nominal  test  conditions  were  an  edge  Mach  number  of  2.9,  nozzle 

stagnation  pressures  from  1  tc  4  atm. ,  and  nozzle  stagnation 

temperatures  in  che  range  of  582  to  592°R  (see  later  discussion 

of  balance  temperature  sensitivity) ,  which  yielded  a  free-stream 

unit  Reynolds  number  range  of  2  to  8  x  106/ft.  For  these 

test  conditions  and  a  constant  wall  temperature  of  532°R,  a 

test  surface  recovery  factor  of  0.86  was  measured  using  a 
thermopile  heat-transfer  sensor;  wall  temperatures  were  held 

uniform  and  constant  by  circulating  a  coolant  (water)  through 
the  test  surface  interior.  All  roughness  tests  were  run  under 
adiabatic  flow  conditions. 
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In  order  to  adapt  the  present  study  to  this  facility 
without  incorporating  major  modifications,  the  apparatus  shown 
schematically  in  Figure  2  was  designed.  In  this  manner,  the 
entire  model  support,  balance,  probe  drive,  etc.,  could  be 
mounted  to  the  original  test  surface  through  existing  instru¬ 
mentation  ports.  The  overall  height  of  this  apparatus,  above 
the  original  test  surface,  was  the  order  of  one  boundary-layer 
thickness  (1.0  inch),  and  it  completely  spanned  the  channel. 

The  model  support  was  located  within  the  constant  pressure 
region  of  the  nozzle.  A  smooth,  shallow  five-degree  ramp, 
followed  by  a  five-degree  expansion  corner  generated  a  new  test 
rhombus  within  which  the  roughness  samples  could  be  mounted. 

The  total  pressure  loss  across  the  weak  oblique  shock  wave  was 
minimal  (the  order  of  1/2  percent  of  nozzle  stagnation  pressure) . 

As  noted  in  Figure  2,  the  full-span  roughness  samples 
extended  over  a  three-foot  axial  length,  beginning  at  the 
expansion  corner.  The  floating  element  portion  of  each  rough¬ 
ness  sample  was  one  foot  in  length,  0.5  foot  in  width  (centered 
about  the  tunnel  centerline)  and  its  forward  edge  was  located 
one  foot  downstream  of  the  expansion  corner.  In  terms  of 
smooth  wall  boundary-layer  thicknesses,  the  floating  element 
covered  an  axial  distance  of  12.4  to  24.8  <5  downstream  of 
roughness  initiation  at  one  atm.,  while  at  four  atm.  the 
corresponding  axial  distances  were  15.7  to  31.4  6,  Profiles 
were  measured  on  centerline,  at  an  axial  station  ~1.5  inches 
upstream  of  the  floating  element  trailing  edge,  corresponding 
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to  a  run  length  over  the  roughness  samples  of  23.2  6  at  one  atm. 
(the  condition  at  which  all  rough  wall  profiles,  except  one, 
were  obtained) . 

Nozzle  stagnation  pressure  and  temperature,  wall  pressure 
and  temperature,  surface  shear  stress,  and  at  least  one  Pitot 
pressure  profile  were  measured  for  each  roughness  model  (shear 
stress  measurements  were  made  in  the  absence  of  any  probes; 
all  profile  data  were  obtained  during  separate  runs) .  Limited 
recovery  temperature  profiles  were  measured  to  verify  certain 
assumptions  used  in  data  reduction  (see  Section  IV, C) .  All 
data  were  monitored  and  recorded  on  the  existing  NOL  Boundary 

/  O  C  \ 

Layer  Channel  data  system v  . 

Figure  3  shows  a  close-up  view  of  the  smooth  wall  floating 
element  section,  with  calibration  weight  applied,  as  installed 
vertically  in  the  tunnel;  the  flattened  tip  (0.005"inch  outside 
thickness)  Pitot  pressure  probe  can  also  be  seen,  as  installed. 
Figure  4  shows  a  close-up  view  of  the  dual  Pitot/recovery 
temperature  probe  used  for  the  simultaneous  measurement  of 

these  variables;  the  temperature  probe  is  based  on  the  work  of 

,,  .  (40) 

Danberg 

Static  pressure  taps  around  the  fixed  periphery  of  the 
smooth  wall  floating  element  were  used  to  check  on  flow 
uniformity.  This  technique  is  by  no  means  conclusive;  however, 
measured  surface  pressure  distributions  indicated  essentially  a 
zero  pressure  gradient  flow,  as  desired,  with  spanwise  Mach 
number  uniformity  within  one  percent  of  nominal  and  axial  Mach 
number  uniformity  within  two  percent  of  nominal.  Centerline  edge 
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Mach  numbers  as  determined  from  measured  wall  static  and  nozzle 
stagnation  pressures  were  consistent  with  those  values  determined 
from  measured  free-stream  Pitot  and  nozzle  stagnation  pressures. 
The  degree  of  flow  two-dimensionality  achievable  within  this 
facility,  for  zero  and  nonzero  pressure  gradient",  has  been 
addressed  in  more  detail  previously  (e.g,,  ref.  (39)). 

A  detailed,  cross-sectional  schematic  of  the  skin-friction 
balance  is  shown  in  Figure  5.  Several  key  points  concerning 
this  instrument,  and  its  accuracy,  bear  mentioning  here. 

2 

The  large  physical  size  of  the  floating  element  (0.5  ft  ) 
was  necessitated  by  the  requirement  to  obtain  integrated  shear 
stress  measurements  over  several  of  the  longest  wavelengths  in 
question. 

Each  floating  element  sample  was  carefully  aligned  with  its 

surroundings  via  the  application  of  selected  shims  to  the 

common  surface  between  sample  and  balance;  final  alignment  was 

verified  by  traversing  a  dial  indicator/stylus  across  the  gap 

between  fixed  and  floating  surfaces  and  by  traversing  a  0.005- 

inch  feeler  gage  around  the  inside  gap  periphery.  In  addition, 

rough/wavy  samples  were  fabricated  and  mounted  such  that  all 

waveforms  possessed  continuity  of  slope  at  the  floating  element 

gap.  Proper  element  alignment  has  been  recognized  as  essential 

(41) 

to  the  generation  of  accurate  skin-friction  data 

This  balance  was  designed  as  a  set  of  parallel  steel  planes, 
joined  together  by  three  pairs  of  steel  webs;  the  entire 
structure  was  fabricated  from  a  single  piece  cf  metal,  i.e.,  the 
two  planes  and  the  webs  were  an  integral  unit.  This  design 
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allowed  for  deflections  in  the  flow  direction  while  maintaining 
the  floating  element  everywhere  parallel  to  the  fixed  baseplate 
(i.e.,  no  pitching  of  the  leading  edge  into  and/or  out  of  the 
flow  was  possible) . 

The  balance  was  mounted  to  the  tunnel  structure  at  a  single 
point,  i.e.,  contact  occurred  only  on  the  annular  surface  between 
the  mounting  shaft  and  mounting  block.  In  addition,  the 
facility  was  designed  to  include  a  plenum  chamber  between  the 
test  surface  and  the  outer  tunnel  wall,  which,  during  testing, 
is  evacuated  to  test  section  pressure.  These  design  features 
isolated  the  balance  from  loads  which  might  have  otherwise  been 
transferred  to  it  from  the  supporting  structure. 

The  center  pair  of  webs  was  instrumented  (both  surfaces) 

with  a  series  of  temperature-compensated  strain  gages;  the 

combined  output  signal  of  these  gages  was  used  to  deduce  surface 

shear  loading.  A  calibration  (load  versus  millivolt  output)  was 

conducted  prior  to  each  run  via  the  application  of  calibration 

weights  (recall  Fig.  3;  a  vertically  oriented  facility  greatly 

aided  in  this  procedure) ,  A  small  hole  was  machined  in  the 

exact  center  of  each  floating  element,  in  which  a  cylindrical 

pin  could  be  mounted;  known  weights  were  then  hung  from  this 

pin.  Excellent  linearity  and  repeatability  were  always  observed 

2 

during  calibration  (loads  to  ''2.5  #^,  shears  to  ~5.0  #^/ft  ). 
After  calibration,  the  pin  was  removed,  its  hole  filled  flush 
with  a  plaster  compound,  and  the  tunnel  sealed  for  running. 
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Calibrations  were  always  conducted  with  the  balance  and 
model  support  at  room  temperature;  no  coolant  passages  existed 
through  th^se  portions  of  the  apparatus.  Initial  tests  were  run 
at  a  nozzle  stagnation  temperature  which,  considering  the  measured 
recovery  factor  of  0.86,  yielded  an  adiabatic  wall  temperature  of 
532°R.  When  initial  wall  (i.e.,  test  sample  surface)  temperature 
differed  by  more  than  several  degrees  from  this  preselected 
adiabatic  wall  temperature,  nonneg] igxble  heating  and/or  cooling 
of  the  test  sample  occurred  prior  to  its  reaching  equilibrium 
(wall  temperature  was  monitored  continuously  before,  during,  and 
immediately  after  each  run;  unfortunately,  room  temperature  could 
not  be  held  constant  over  a  long  time  scale,  i.e.,  from  one 
run  to  the  next) .  For  such  cases,  unacceptable  zero  shifts 
were  noted  in  the  balance  output  signal  (here  defined  as  >2% 
of  full-scale  load) .  A  test  procedure  was  then  formulated  to 
circumvent  this  problem;  nozzle  stagnation  temperature  was  pre¬ 
selected  and  accurately  maintained  (within  ~2°R,  via  a  feedback 
control  system  on  the  flow  heater)  such  that  adiabatic  wall 
temperature,  for  each  particular  run,  equalled  the  prerun  wall 
temperature.  In  this  manner,  no  appreciable  thermal  effects  were 
imposed  on  the  balance  and  observed  zero  shifts  were  thus  held 
within  acceptable  limits  (<2%  of  full-scale  load;  actually, 
most  observations  were  within  1%).  This,  coupled  with  the 
levels  of  linearity  and  repeatability  achieved  during  calibration, 
resulted  in  the  following  estimates  for  shear  load  accuracy: 

At  the  highest  R8q  levels  (4  atm.),  where  the  highest  levels  of 
shear  were  experienced,  measurement  accuracy  was  of  the  order  of 
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three  percent;  at  the  lower  ReQ  levels  (1  atm.),  where  reduced 
shear  levels  were  experienced,  data  accuracy  was  of  the  order 
of  ten  percent. 

IV  RESULTS  AND  COMPARISONS  WITH  EXISTING  DATA 
A.  Smooth  Wall  Boundary-Layer  Characteristics 

Before  presenting  and  discussing  the  bulk  of  the  skin- 
friction  measurements,  characteristic  dimensions  of  the  reference, 
or  smooth  wall,  boundary  layer  should  be  noted.  Figure  6  shows 
momentum,  mass  displacement,  and  boundary-layer  thickness  as  a 
function  of  free-stream  unit  Reynolds  number.  For  comparison 
purposes,  the  minimum  measured  smooth  wall  momentum  thickness 
was  still  nearly  double  the  maximum  roughness  dimension  and/or 
maximum  wave  amplitude  tested. 

The  smooth  wall  boundary- layer  scale  of  real  importance 
for  any  study  of  roughness  effects  is,  however,  the  laminar,  or 
viscous,  sublayer  thickness.  By  combining  the  definition  of 


friction  velocity. 


»  JZ 

T  TP 

*  r  TaT 


with  Newton's  law  of  friction, 


or,  for  small  y, 


,3u. 
T  "  11  <Sy> 


Tw  °  >J„lu/y> 


the  sublayer  velocity  distribution  can  be  shown  to  be  stated  by 


v  ■  (?) 


+  + 
u  =  y 
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Based  on  experimental  evidence,  the  outer  edge  of  the 
viscous  sublayer  occurs  at  a  y+  value  -11.0.  Thus, 


11  v 


5s  ”  u. 


Introducing  the  definitions  of  v  ,  u  ,  ,  c  and  M  into  Equation 

(6) ,  coupled  with  the  assumptions  of  ideal  gas  and  p^  =  p 


results  in 


11  p  c 

6  =  - Z-Z-  (7) 

s 

Thus,  6  can  be  determined,  for  known  P  ,  T  values  from 
s  c  '  oM 

00  CO 

surface  measurements  of  p  ,  T  and  x  .  Results  so  generated 

*w'  w  w  3 

are  shown  in  Figure  7. 

The  necessity  for  determining  viscous  sublayer  thicknesses 

from  surface  measurements  is  best  illustrated  when  one  views 

an  overlay  of  the  present  probe-tip  dimensions  on  Figure  7 . 

For  an  overall  tip  dimension  of  0.005  inch,  the  closest  point 

to  the  wall  at  which  a  velocity  measurement  can  be  claimed  is 

0,0025  inch,  or  at  one-half  the  probe-tip  height.  The  problem 

of  probe  tip/wall  interference  must  then  be  addressed. 

Figures  8a-d  show  present  smooth  wall  velocity  profiles  in 

terms  of  (u/u^) ,  (y/6)  coordinates.  Both  uncorrected  and 

corrected  data  points  are  shown  (probe  corrections  were  based  on 

(39) 

earlier  NOL  research,  e.g.,  Voisinet  ,  and  were  applied  to 
present  profile  data  only  where  applicable,  i.e.,  only  for  smooth 
wall  conditions) .  As  can  be  seen,  these  probe  corrections 
influenced  the  data  only  in  the  immediate  vicinity  of  the  wall; 
integral  parameters  shown  earlier  in  Figure  6  were  computed  from 
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the  corrected  velocity  distributions,  but  any  differences 
between  those  values  shown  in  Figure  6  and  values  computed  from 
the  uncorrected  velocity  distributions  were  indistinguishable 
(third  significant  figure  influence  only) . 

One  way  to  demonstrate  self-consistency  between  present 
surface  shear  and  velocity  profile  data  is  to  plot  sublayer 
velocity  distributions,  as  determined  from  measured  shears,  on 


Figures  8a-d.  From  Equation  (3) , 


y  u  "I 

T  = 

w  L  6  J 


tu/uj 

(y/fi) 


(VuJ 


[V.1 


Cy/fi) 


On  a  log-log  plot  in  these  coordinates,  the  sublayer  velocity 
distribution  will  appear  as  a  straight  line  of  slope  equal  to 

45  degrees,  i.e.,  from  Equation  (9), 

T  t  6  ~\ 


log1()(u/uj  «  log. 


r  x  6  i 

w 

P  U 

L  W  ooj 


+  log10  (y/,!) 


u'  =  constant  +  y‘ 


1  r45°  slope] 


and  its  intercept  with  the  vertical  axis  used  on  Figures  8  will 


occur  at 


intercept  ^ 


•°oiifey 


Results  so  plotted  show  that,  in  every  case,  the  probe- 
corrected  velocity  distribution  approached  the  shear-determined 
sublayer  velocity  distribution  in  the  vicinity  of  the  wall. 
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Also,  in  every  case,  an  extrapolation  of  the  logarithmic  region 

slope  (profile  data]  to  its  intercept  with  the  sublayer  velocity 

distribution  (shear  data]  defined  a  point  consistent  with  the 
edge  of  the  viscous  sublayer,  as  determined  from  measured 
shear. 

Smooth  wall  temperature  and  velocity  profile  data  will  be 
discussed  in  more  detail  in  Sections  IV,  C  and  D,  but  first  a 
presentation  of  all  smooth  and  rough  wall  shear  stress  data 
will  be  made. 

B.  Surface  Shear  Stress  Measurements 

Figures  9a-c  summarize  the  present  skin-friction  coefficient 
results  as  a  function  of  smooth  wall  Re0 .  The  data  are  sub¬ 
divided  into  these  groupings  for  ease  in  interpretation.  Smooth 
wall  coefficients  are  plotted  on  each  graph  for  reference 

purposes.  The  smooth  wall  skin-friction  coefficient  was  found 

-0  14 

to  decrease  with  increasing  Re0,  such  that  *  Re0 

o 

Exponents  in  the  range  of  -0.10  to  -0,20  are  expected  for  fully 
developed  turbulent  boundary  layers.  Repeatability  of  these 
measurements  is  illustrated  by  those  two  cases  (smooth  and 
molded  models)  wherein  the  models  were  remounted,  realigned, 
and  tested  a  second  time. 

Figure  9a  shows  the  sand -grain  results.  These  data  are 
quite  similar  in  functional  dependence  to  the  Mach  3  sand-grain 
results  of  Goddard  (see  Fig.  18  of  ref.  (9)).  The  80-grit 
data  serve  especially  well  to  illustrate  the  relative  importance 
of  the  two  dominant  physical  scales  involved,  roughness  height 
and  viscous  sublayer  thickness.  At  one  atmosphere  total  pressure, 
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the  viscous  sublayer  thickness  was  seven  mils  {Fig.  7)  ,  while 
roughness  height  was  six  mils,  i.e.,  roughness  peaks  had  not 
yet  begun  to  protrude  outside  the  sublayer.  Correspondingly, 
no  increase  in  skin-friction  coefficient  above  the  measured 
smooth  wall  value  was  observed.  However,  as  Reynolds  number 
was  increased  the  viscous  sublayer  thinned  while  roughness 
dimension  remained  constant.  The  net  result  was  an  increase  in 
the  rough  wall  skin-friction  coefficient  above  the  smooth  wall 
level.  Roughness  effects  are  thus  generally  subdivided  into 
three  regimes: 

for  k  <  6  ,  Re.  <  10  or  11;  aerodynamical ly  smooth 

for  <  k  <  6  or  7  {  ,  11  <  Re,  <  '-70;  transitionally  rough 

S  <3  K 

for  k  >  6  or  7  6  ,  Re.  >  -70  ;  fully  rough 

S  K 

where  ku 

Re^  =  =  roughness  Reynolds  number  (12) 

Figure  9b  presents  a  subset  of  data  which  shows  effects  of 
superimposing  roughness,  separately,  on  a  short,  and  on  a  long 
periodic  waveform,  as  well  as  the  effects  of  superimposing 
roughness  on  a  surface  which  simultaneously  possesses  both  a 
short  and  a  long  periodic  waveform.  For  these  models,  k,  egw, 

^SW'  eLW  anc^  ^LW  were  held  constant. 

Roughness  as  well  as  roughness  superimposed  on  a  short 
periodic  waveform  yielded  essentially  the  same  result.  Rough¬ 
ness  superimposed  on  a  long  periodic  waveform  actually  showed 
some  reduction  in  surface  shear  as  compared  to  roughness  alone. 
The  important  point,  however,  is  that  superposition  of  a  given 
roughness  dimension  on  a  surface  which  simultaneously  possessed 
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both  short  and  long  periodic  waveforms  increased  surface  shear 
above  .levels  measured  for  identically  rough  surfaces  possessing 
only  one,  or  "•''ne,  of  these  same  periodic  waveforms. 

Figure  9c  presents  the  final  data  subset,  consisting  of  the 
motor  case  sample  and  three  machined  samples  which  simultaneously 
possessed  roughness,  as  well  as  short  and  long  wavelengths.  It 
was  hoped  that  these  three  machined  samples  would  simulate  the 
motor  case  material  over  different  portions  of  its  boost 
trajectory  (by  varying  k,  e  and  \  scales  versus  sublayer 
thickness) ,  As  can  be  seen,  the  desired  simulations  were  not 
achieved  (models  9  and  10  should  have  given  nearly  identical 
results  had  accurate  simulation  been  achieved) .  This  could  be 
due,  in  part,  to  the  fact  that  the  motor  case  material  was  of 
a  random  rough/wavy  character  (a  fact  which  should  be  kept  in 
mind,  even  though  analysis  of  surface  contour  traces  resulted 
in  definition  of  prominent,  or  reoccurring,  roughness  and 
waveform  dimensions,  i.e.,  those  listed  in  Table  1). 

Increasing  roughness  and  waveform  dimensions,  while 
maintaining  nearly  constant  (e/A)  ratios  (-0.01  for  short 
waves  and  ~ 0.005  for  long  waves)  resulted  in  increased  surface 
shear  stresses.  All  three  samples  showed  a  similar  functional 
dependence  on  ReQ.  On  the  other  hand,  the  motor  case  material 
showed  an  opposite  trend,  with  decreasing  slightly  as  Re0 
increased. 

Equivalent  sand-grain  roughnesses  for  the  motor  case 
material  were  obtained  by  cross -plotting  rough- to -smooth 
wall  ratios  versus  sand-grain  height,  and  superimposing 
motor  case  ratios  on  corresponding  constant  Reynolds  number 
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(i.e.,  constant  P  )  lines.  The  results  are  shown  in  Figure  10. 

U(X) 

Equivalent  sand-grain  roughnesses  for  this  material  (5  to  10 
mils)  were  3.8  to  7.6  times  t-ne  roughness  scale  attributed  to 
it,  and  1.25  to  2.5  times  its  ascribed  short  wave  amplitude. 

As  suspected,  skin-friction  drag  calculations  based  solely  on 
the  smallest  roughness  scale  were  not  representative  of  the 
actual  situation. 

This  section  is  concluded  with  two  comparisons  between 
present  results  and  previously  published  data.  Figure  11  shows 
the  first  comparison,  strictly  for  smooth  wall  data.  As  can  be 
seen,  present  smooth  vail  skin-friction  coefficients  are  in 

reasonable  agreement  with  the  findings  of  previous  investigators 

(42) 

and  with  the  correlation/calculation  scheme  of  Spalding-Chi 
(within  "20%  of  the  latter,  as  calculated  for  a  nominal  Mach 

(43) 

'lumber  of  3;  data  and  calculations  are  as  reported  by  Sturek  ). 
Some  recent,  unpublished  data  of  Voisinet,  obtained  on  the  flat 
nozzle  wall  of  the  NOL  Boundary  Layer  Channel,  via  a  small 
floating  element  balance,  are  also  shown. 

A  second,  more  pertinent,  comparison  is  shown  in  Figure  12, 
wherein  the  rough-to-smooth  wall  skin-friction  coefficient  ratio 
is  plotted  as  a  function  of  the  logarithm  of  roughness  Reynolds 
number.  These  correlating  parameters  were  originally  proposed 
by  Nikuradse v  1  (incompressible  flow)  and  Goddard  (compressible 

flow)  and  were  based  on  their  sand-grain  results.  The  present 
data  plus  those  Y-groove  roughness  results  of  Young^^,  Mann^°^, 
and  Wade^1^,  are  shown  in  comparison  with  the  bounds  of 
Goddard's  ^  compressible,  sand-gr  in  results  (note  that  all 
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data  shown  are  for  adiabatic  flow) ,  Several  key  points  are 
mentioned  in  the  succeeding  paragraphs. 

(9) 

1.  Present  sand-grain  results  resubstantiate  Goddard's 
correlation  and  the  conclusion  emanating  from  it,  namely  that 
’'the  effect  of  surface  roughness  on  skin-friction  drag  is 
localized  deep  within  the  boundary  layer  at  the  surface  itself 

and  is  independent  of  the  external  flow,  i.e.,  Mach  number,  per  se, 

is  eliminated  as  a  variable."  The  significance  of  the  breaking 

point  in  the  correlation,  (C^/C^  )  >  1  for  Re^  >  10  or  11,  has 

o 

been  discussed  earlier  in  this  section  (Eqs.  (12)). 

2.  A  similar  functional  dependence  between  (Cf/C^  ) 

o 

and  Re^  has  also  been  demonstrated  for  machined  (V-groove) 
roughnesses,  as  superimposed  on  both  flat  and  single  periodic 
waveform  surfaces,  via  present  data  and  previously  published  results. 
For  these  type  surfaces,  equivalent  sand-grain  roughnesses  must 
not  differ  substantially  from  actual  roughness  dimensions, 
considering  the  observed  level  of  agreement  with  the  bounds  of 
Goddard's  sand-g._ain  results, 

3.  Superposition  of  roughness  on  surfaces  which 

simultaneously  possessed  both  short  and  long  periodic  waveforms 

(models  10,  11,  12)  resulted  in  a  narrow,  self-consistent  band 

of  data,  but  one  which  failed  to  fall  within  the  bounds  of 

Goddard's  correlation.  Rather,  these  results  (plotted  using  k, 

not  £gW/  as  the  roughness  scale  in  Re^)  ,  departed  from  the 

(Cf/Cf  )  =1.0  line  at  Re^  values  near  4  to  5 ,  and  thereafter 
o 

exhibited  a  functional  dependence  on  Re^  similar  to  that  of  the 
original  correlation.  Such  behavior  is  indicative  of  the 
amplification  in  surface  shear  stress  caused  by  this  type  of 
rough/multiple  waveform  surface.  Increases  in  surface  shear. 
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above  smooth  wall  levels,  were  observed  for  these  cases  even 

though  the  superimposed  roughness  scale,  k,  would  not  have  been 

sufficient,  by  itself,  to  cause  such  increases,  i.e.,  at  the 

departure  from  the  (C^/C^  )  =  1.0  line  the  roughness  dimension  of 

o 

each  rough/multiple  waveform  surface  was  always  less  than  one- 

half  the  corresponding  sublayer  thickness. 

4.  The  motor  case  (i.e.,  random  rough/wavy  wall)  data, 

plotted  using  its  ascribed  k  value,  did  not  correlate  with  any  of 

the  other  results.  A  horizontal  shift  applied  to  these  data  via 

a  replacement  of  the  k  roughness  scale  in  Re^  with  e ,w  (ascribed 

short-wave  amplitude)  brought  them  within  the  bounds  of  other 

results,  but  the  weaker  functional  dependence  on  roughness  Reynolds 

number  remained,  of  course,  unchanged,  i.e.,  the  slope  of  (C^/Cf  ) 

o 

versus  Re^  remained  below  that  possessed  by  other  results. 

C.  Smooth  and  Rough  Wall  Temperature  Profiles 

Limited  temperature  profile  data  were  obtained  for  smooth 

wall  and  #24  grit  rough  wall  conditions.  These  data  were  needed 

to  verify  certain  temperature-velocity  relationships  frequently 

assumed  in  compressible,  adiabatic  flow  data  reduction.  This 

section  summarizes  those  results . 

The  generalized  total  temperature-velocity  relationship 
(44) 

of  Danberg  can  be  written, 

T  =  (3u  +  Cl  -  (3 ) u 2  (13) 


where , 


(recovery  factor) 
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This  generalised  expression  reduces  to  two  well-known 
special  cases.  First,  for  unit  Prandtl  number,  where  r  =  1.0 


T  =  T 
aw  o. 


B  =  1.0 


and, 

T  =  u  (Crocco  relation) 
Second,  for  adiabatic  flow,  where'r  /  1.0*, 


(14) 


T  =  T 
w  aw 


6  =  0 


and, 


—  —2 

T  =  u  (quadratic  or  Walz  relation) 


(15) 


Figure  13  presents  smooth  and  #24  grit  rough  wall  data  in 
comparison  with  Equations  (14)  and  (15).  A  total  temperature 
overshoot,  and  agreement  with  the  quadratic  expression,  were 
noted  in  each  case,  characteristic  of  adiabatic,  nozzle  wall 
boundary  layers. 

For  use  in  actual  profile  data  reduction,  a  static 
temperature-velocity  relationship  was  desired.  According  to 
Walz v  ,  a  generalized  static  temperature-velocity  relationship 
can  be  written  as, 


(±H  =  A  +  B  (H-l  +  C  (£-} 


,u 


u 


u 


(16) 


where  coefficients  are  evaluated  from  the  following  boundary 
conditions : 

(1)  for  u  =  0  ,  T  «  Tw 

(2)  for  u  =  uTO  ,  T  =  Tro 

3  T 

(3)  for  adiabatic  flow,  T  =  T  „  and  (s~)  =  0;  thus, 

v  av  dY  w 


*  Simultaneous  assumption  of  r  =  1,0  and  adiabatic  flow  would 

leave  (5  undefined,  i.e.,  T  .  =  T  =  T  and  0  =  0/0. 

aw  o„  w 
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,3T, 

C5y’, 


0 


w 


but. 


l3yJ 


v 


?  0, 


therefore , 


(2Zj 


0 


where  u  =  0 


Thus : 


(r  *  1.0) 


For  adiabatic  flow,  where  6  «=  0  *=  B,  Equations  (13)  and  (16) 
can  be  shown,  mathematically,  to  be  equivalent.  Figures  14a-b  show 
present  results  in  comparison  with  Equation  (16) .  As  can  be 
seen,  the  agreement  is  very  good.  Consequently,  Equation  (16) 
was  used  to  reduce  all  profile  data  where  actual  temperature 
measurements  were  not  made. 

D,  Smooth  and  Rough  Wall  Velocity  Profiles 

In  order  to  effectively  present  and  discuss  velocity  profile 
results,  certain  mathematical  considerations  for  smooth  and  rough 
wall  turbulent  boundary  layers  must  be  reviewed. 

Since  we  are  dealing  with  compressible  flow,  a  velocity 
transformation  to  the  incompressible  plane  is  required  before  any 
comparisons  between  present  results  and  previously  published 
mathematical  formulations  and/or  experimental  results  can  be 
undertaken.  In  this  investigation,  a  transformation  based  on 
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the  work  of  Van  Driest^*  '  was  used': 


*  -  (U<°\ 
u  "  U  / 


+  sm 


where , 


(b2  +  4A2 


T  -  T 
aw  ® 


•{— -s— ' ) 

\Vb2  +  4A2/ 


’  T 

aw  . 

T  1 

L  w  J 


r  =  recovery  factor  ^  1.0 

Velocities  appearing  in  all  subsequent  correlation  schemes 
(i.e.,  [u+,y+],  t  (u  -  u^)/^,  (y/S)]  and  [(Au/ut),  Rekl  coordinates) 
are  transformed  values;  continued  use  of  the  *  superscript  will 
not  be  made. 

An  excellent  review  of  smooth  and  rough  wall  turbulent 

(3) 

boundary-layer  correlation  techniques  was  given  by  Clauser  ; 
a  brief  summary  follows. 

Within  the  turbulent  boundary  layer  three  distinct  regions 
have  been  found  to  exist: 

1.  an  inner,  or  viscous  sublayer  region,  where 


...  • .(?) 

'  w  > 


2.  an  outer,  or  velocity  defect  region  (sometimes 
referred  to  as  the  wake  region) ,  where 


’u  ~  U»1 
-  UT  J 


f  (y/<5) 


*  As  programmed  in  the  NOL  Law  of  the  Wall  -  Law  of  the  Wake  data- 
reduction  program;  the  assistance  of  Robert  L.  P.  Voisinet  in 
adapting  this  program  for  use  in  the  present  investigation  is 


gratefully  acknowledged. 
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3.  an  overlapping,  or  law-of-the-wall  region. 

In  order  to  deduce  a  mathematical  formulation  for  the  over¬ 
lapping  region  (Eq,  (4)  describes  the  sublayer) ,  the  following 
argument  was  presented:  Within  the  overlapping  region,  both 
Equations  (18)  and  (19)  must  hold;  therefore,  equating  (u/u^) 


expressions , 


<y/‘>  ♦  (£)  -  4  ■  P) 


A  comparison  of  the  two  sides  of  Equation  (20)  shows  that  the 

/«M 

effect,  of  the  multiplication  factor  I- — 1  inside  g  must  be 


equivalent  to  the  additive  term 


ft) 


inside  g  must  be 


outside  f;  the  logarithm 


is  the  only  function  with  this  property.  Thus,  within  the  over¬ 
lapping  region. 


(u  -  uJ/uT  -  A  log1Q(y/6)  +  B  (21) 

(u/uT)  =  A  lOg10^J^  +  C 

u+  =  A  log10(y+)  +  C 

Equation  (22'  defines  the  law-of-the-wall ,  in  the  absence  of 
roughness  effects.  Equality  of  slopes  for  Equations  (21)  and 
(22),  within  their  mutual  regime  of  validity,  can  be  demonstrated 
by  equating  (u/uT)  expressions  from  each  and  differentiating 
both  sides  with  respect  to  the  nondimensional  coordinate  (y/6). 
This  requirement  for  equalicy  of  slopes  within  the  law-of-the- 
wall  region  provides  one  self-consistency  check  on  profile  data. 

Coles  has  formulated  an  additive  term  for  Equation  (22) 
in  order  to  generate  a  simultaneous  mathematical  expression  for 
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both  the  law-of-the-wall  and  law-of -the -wake  regions  (i.e.  .  for 
all  y  values  >  11).  However,  since  primary  emphasis  here  has 
been  placed  on  roughness  effects,  and  since  roughness  does  not 
strongly  influence  the  outer  or  wake  region  (as  will  be  demon¬ 
strated  in  Fig.  17),  no  attempts  have  been  made  to  compare 

present  profile  data  with  Coles’  formulation  (previous  NOL 
(39' 

studies  ' ,  involving  pressure  gradient  effects  on  smooth  wall 
turbulent  boundary- layer  development,  have  included  such  analyses) . 

Roughness  effects  on  the  mean  velocity  profile  are  primarily 
concentrated  in  the  inner  20  percent  of  the  layer,  i.e.,  within 
the  viscous  sublayer  and  law-of-the-wall  (logarithmic)  region. 
Effects  of  roughness  on  the  law-of-the-wall  have  been  reasoned, 
and  verified  experimentally,  to  be  reflected  solely  as  a  shift 
in  the  intercept  C  of  Equation  (22) .  The  shift,  itself,  is  a 
function  of  the  roughness  Reynolds  number,  Rek,  i.e.,  for  rough 
walls , 

<u/uT)  =  A  lCio(y  +  c  -  (^)  (23) 

where , 

+  + 

A  =  smooth  wall  u  ,y  slope 

(57)  =  f (Rek>  (24) 

For  large  roughness  Reynolds  numbers  (Re^  >  ~70)  ,  the  sub¬ 
layer  is  destroyed  and  the  inner  regions  of  the  boundary  layer 
must  become  independent  of  viscosity.  For  these  conditions  to  be 
met,  Equation  (24)  must  be  of  the  form 

/ku  \ 

=  A  log1Q  — 1  +  D  [fully  rough]  (25) 
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which  can  bo  verified  by  substituting  Equation  (25)  into 
Equation  (23)  -■  the  resulting,  fully  rough,  expression 

u+  =  A  log10(y/k)  +  C  -  D  Ifully  rough] 


shows  no  viscosity  dependence. 


(26) 


The  fact  that  the  slope  of  CAu/ut)  versus  log^Q(Ret),  for 


Re^  >  ~70,  mist  match  the  smooth  wall  u  ,y  slope,  provides 
another  self-consistency  check  on  profile  data. 

A  technique  originated  by  Clausar^^  provides  still  another 
self-consistency  check  between  measured  shear  and  profile  data. 
Here  the  surface  shear  stress  is  inferred  from  the  velocity 
profile  slope.  Eliminating  t  from  the  definitions  of  uT  and 


results  in 


UT  = 


V 


Cf  Poo 


W 


U 


(27) 


Substituting  Equation  C27)  into  Equation  (23)  and  rearranging 
yie Ids, 

V~Cf  pm  ' 

2 - —  log^pCy)  +  tterms  with  no  y  dependence]  (28) 

"  pw 

This  technique  then  calls  for  plotting  (u/u^)  versus  log^^ty) 

and  measuring  the  slope,  S,  in  the  logarithmic  region.  From 
Equation  (28) 

S  =  A^l  ~  ~~  (29) 


or, 


T1 

Koo 

!  2 

pw 

s2 

pw 

A2 

Pro 

(30) 
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Inferred  skin-friction  coefficients  were  obtained  in  this 
manner  for  two  limiting  cases  and  the  results  are  shown  below: 


meas. 


^infer . 


'meas. 


smooth  wall, 

Cl  atm) 

fully  rough  wall, 
(#24  grit,  2  atm.) 


1.81  x  10  1.69  x  10 


3.51  x  10  3.17  x  10 


-0.066 


-0.097 


Levels  of  agreement,  within  measurement  accuracy,  were  noted. 

In  order  to  conduct  any  such  analysis  of  rough  wall  profile 
data,  the  question  of  effective  y  origin  must  be  addressed. 

Figure  15  schematically  depicts  this  problem.  Actual  Pitot 
probe  measurements  can  only  be  made  down  to  the  roughness  peaks , 
i.e.,  to  the  pcint  of  probe  contact.  However,  previous  investi¬ 
gators  (e.g.,  Srottron^13^  and  Perry ^14^)  have  found  that  the 
effective  y  origin  lies  somewhere  between  the  roughness  peaks 
and  valleys . 

Systematic  variations  were  therefore  applied  to  the  y  origin, 
such  that 


where , 


y'  (shifted  origin)  =  y  (probe  origin)  +  ak 


a  =  0,  1/2,  1 


v  -  V  e  *  p  * 

'  SW'  'LW 


+  +. 


and  corresponding  velocity  profiles  were  plotted  in  both  [u  ,y  ] 
and  I  (u  -  uro)/uT,  (y/<5)]  coordinates,  e.g.,  Figure  16.  An 
effective  y  origin  was  thereby  defined  by  that  particular  shift 
(a  •  kj  which  brought  the  rough  wall  velocity  profile  data  parallel 


*Models  7  through  12,  where  appropriate. 
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+  + 

to  the  smooth  wall  u  ,y  slope,  within  the  logarithmic  region 
(recall  Eq.  (23)).  Examination  of  all  such  plots  showed  that 
the  combination  of  a  =  1/2,  k  =  k,  most  consistently  met  this 
criterion,  even  for  the  rough/wavy  wall  cases  (evidently,  those 
larger  shifts  associated  with  waveform  amplitudes  were  un¬ 
realistic)  .  The  #24  grit,  1  atm.,  profile  data  were  selected  for 
presentation  in  Figure  16  because  they  most  clearly  illustrate 
the  features  discussed  above. 

An  effective  y  crigin,  determined  in  this  manner,  is  utilized 

in  all  subsequent  profile  plots.  Note  that  the  continued  use 

of  the  prime  superscript  on  y  will  not  be  made. 

Fiqure  17  shows  present  profile  data,  in  terms  of  velocity 

defect  coordinates,  in  comparison  with  the  range  of  smooth  and 

(3) 

rough  wall  incompressible  results  reported  by  Clauser  .  The 
level  of  agreement  is  seen  to  be  quite  good.  Several  points 
concerning  this  figure  should  be  made. 

1.  All  smooth  and  rough  wall  data  collapse  to  a 
nearly  universal  cur/e  in  the  outer  80  percent  of  the  layer 
(i.e,,  in  the  wake  region,  (y/6)  >  ~0.2),  showing  that  the  effects 
cf  roughness  are  indeed  localized  deep  within  the  boundary  layer. 

2.  As  one  views  data  below  the  wake  region  (inner 

20  percent  of  the  layer) ,  some  variation  between  smooth  and  rough 
wall  resuits  becomes  apparent.  With  the  exception  of  the 
random  rough/wavy  (molded  surface)  data,  all  rough  wall  data  are 
displaced  slightly  downward  (in  these  coordinates)  from  the 
smooth  wall  curve.  Similar  observations  have  been  made 
elsewhere  . 


27 


V 


NOLTR  74-34 

3.  Ail  smooth  and  rough  wall  data,  within  the 
logarithmic  region,  possess  slopes  which  match  the  smooth  wall 
u  ,y  slope,  thereby  demonstrating  a  self-consistency  check 
noted  earlier. 

Figures  18a-c  show  present  profile  data  plotted  in  terms 
of  u+,y+  coordinates.  Data  are  subdivided  into  three  groups,  as 
was  done  in  Figure  9,  for  ease  in  interpretation.  As  before, 
the  smooth  wall  results  are  shown  on  each  plot  for  reference 
purposes. 

Figure  18a  shows  smooth  and  sand-grain  results.  These  data 
cover  the  entire  range  of  possible  roughness  regimes,  from 
aerodynamically  smooth  to  fully  rough.  The  curved  solid  line 
represents  the  theoretical  sublayer  velocity  distribution 
(u+  =  y+) .  The  straight  solid  line  through  the  smooth  wall  data 
is  described  by  Equation  (22),  i.e,,  the  law-of-the-wall  cor¬ 
relation.  Present  smooth  wall  slope  and  intercept  were  found 
to  be  4.33  and  5.50  (at  1  atm.),  respectively.  Under  the  in¬ 
fluence  of  roughness,  velocity  profiles  were  seen  to  shift  down¬ 
ward,  to  the  right,  from  the  sii  ooth  wall  curve,  while  remaining 
parallel  to  it  within  the  logarithmic  region  (all  straight  solid 
lines  are  parallel) .  Such  observations  resubstantiate  the 
discussions  centered  around  Equation  (23) . 

Less  notable  velocity  shifts  (if  any)  were  observed  for  the 
molded  and  machined  surfaces,  due  primarily  to  conditions  at 
which  measurements  were  made  (i.e.,  at  1  atm.  total  pressure, 
where  most  values  for  these  models  were  found  to  approximate 
the  smooth  wall  level)  ,  However,  results  shown  in  Figures  18b-c 
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lend  themselves  to  an  important  conclusion,  namely,  that  velocity 
profiles  over  rough/wavy  surfaces  (including  random  rough/wavy 
and  rough/multiple  periodic  waveform  surfaces)  possess  logarithmic 
regions  wherein  the  law-of-the-wall  correlation  is  valid. 

Figure  19  concludes  this  section  with  a  comparison  between 
present  and  previously  published  data,  showing  roughness-induced 
velocity  shift  (Au/u  )  as  a  function  of  Re,  .  Observed  velocity 
shifts  for  flows  over  sand-grain  roughened  surfaces  were  found 
to  be  in  close  agreement  with  other  compressible  and  incompressible 
sand-grain  results.  Insufficient  velocity  shift  data  were  obtained 
for  flows  over  other  surface  conditions  to  warrant  any  meaningful 
comparisons  and/or  statements. 

V  CONCLUSIONS 

Based  on  present  results,  their  demonstrated  levels  of  self- 
consistency,  and  favorable  comparisons  with  previously  published 
results,  the  following  conclusions  are  stated: 

1.  Superposition  of  a  given  roughness  dimension  on  a 
surface  which  simultaneously  possessed  both  short  and  long 
(shallow)  periodic  waveforms  increased  surface  shear  above  levels 
measured  for  identically  rough  surfaces  possessing  only  one,  or 
none,  of  these  same  periodic  waveforms. 

2.  Equivalent  sand-grain  roughnesses  for  the  random 
rough/wavy  wall  (motor  case  material)  were  of  the  order  of  four 
to  eight  times  the  physical  roughness  scale  attributed  to  it. 

3.  Goddard's  compressible  skin-friction  correlation, 

(Cf/Cf  )  =  constant  Ilog^Q  Re^l  +  constant,  for  Re^  >  ~  10,  has 
o 

been  resubstantiated  for  sand-grain  roughnesses.  A  similar 
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"V 


functional  dependence  between  these  variables  has  also  been 

demonstrated  for  machined  (V-groove)  roughness  patterns/  as 

applied  to  both  flat  and  single  periodic  waveform  surfaces. 

Skin-friction  data  obtained  on  rough  surfaces  possessing  multiple 

periodic  waveforms  were  self-consistent  and  exhibited  a  like- 

functional  dependence  on  Re^,  but  showed  departure  from  the 

(C-/C-  )  -  1.0  line  at  lower  Re,  values  (~  5  as  opposed  to  -  10 
o  K 

for  all  other  results) . 

4.  Walz'  temperature-velocity  relationship,  for 
adiabatic  flow,  was  found  to  accurately  describe  measured 
temperature  profiles  over  both  smooth  and  rough  surfaces. 

5.  Velocity  profiles  measured  over  rough/wavy  surfaces, 
including  random  rough/wavy  and  rough/multiple  periodic  waveform 
surfaces,  were  found  to  possess  logarithmic  regions  wherein  the 
law-of- the -wall  correlation  was  valid. 

6.  When  plotted  in  terms  of  velocity  defect  coordinates, 
present  smooth,  rough,  and  rough/wavy  profile  data  collapsed  to 

a  near  universal  curve  in  the  outer  portions  of  the  layer,  in 
agreement  with  previously  published  results,-  verifying  that, 
even  for  the  complex  surface  patterns  considered  here,  roughness 
effects  are  localized  deep  within  the  boundary  layer. 
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3  PHOTOGRAPH  OF  FLOATING  ELEMENT,  WITH  CALIBRATION 
WEIGHT  APPLIED,  AND  PITOT  PROBE  INSTALLED 
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FIG.  5  SKIN-FRICTION  BALANCE  SCHEMATIC 
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FIG.  7  SMOOTH  WALL  SUBLAYER  THICKNESSES  FROM  MEASURED  SHEAR 
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FIG.  8,a-b  SMOOTH  WALL  VELOCITY  PROFILES 
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FIG.  9-a-b-c  SMOOTH  AND  ROUGH  WALL  SKIN  FRICTION  COEFFICIENTS 
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SMOOTH  AND  ROUGH  WALL  TEMPERATURE-VELOCITY  RELATIONSHIPS 
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SMOOTH  AND  ROUGH  WALL  TEMPERATURE-VELOCITY  RELATIONSHIPS 
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